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Abstract. The purpose of this work is to construct asymptotically flat, time 
symmetric initial data with an apparent horizon of prescribed intrinsic geome- 
try. To do this, we use the parabolic partial differential equation for prescribing 
scalar curvature. In this equation the horizon geometry is contained within 
the freely specifiable part of the metric. This contrasts with the conformal 
method in which the geometry of the horizon can only be specified up to a 
conformal factor. 



1. Introduction 

In this work, the excision approach to constructing black hole initial data is used. 
Specifically, black hole initial data for the Einstein equations is taken to consist of 
the following: (1) a manifold M = R^\if , where if is a union of a finite number of 
balls. (2) a metric g an. M asymptotic to the flat background metric 5; in this work 
we take this in the very strong sense that there exist constants Cj for all J > such 
that \r^^^dj{gki — 5ki)\ < Cj, where r is the coordinate radius. (3) a second rank 
tensor field k asymptotic to 0. (4) a vectorfield J asymptotic to 0. (5) a scalar 
field p asymptotic to 0. Furthermore, we require that {g, J, fc, p) satisfy the Einstein 
constraint equations 

R{g) + {ivkf - \k\l = 16^p 
V • fc - V(tr/c)2 = -SttJ, 

the dominant energy condition 

P> \ J\g. 

and, in addition, an apparent horizon boundary condition on dK. To write the 
latter, let n, H be the outward unit normal and mean curvature of dK, respectively. 
The apparent horizon boundary condition is that on dK, one has 

0+ = H + k{n, n) - trA: = 

= -H + k{n, n) - trfc < 0. 

In addition, there should be no other compact surfaces satisfying 6*+ < 0. 

The interpretation of the data is that M is a spacelike slice with induced metric 
g and second fundamental form k in an asymptotically flat spacetime. The fields 
J and p are, respectively, the local momentum current density and mass energy 
density according to an observer comoving with M. The conditions that we have 
assumed on dK imply, assuming weak or strong cosmic censorship and the weak or 
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strong energy condition, that if the data were extended inside of dK, this surface 
would form the boundary of the totally trapped region (for a definition see [17j). 
and thus dK is the apparent horizon. The surface dK represents a black hole 
within M since, again, if {M, g, k, J, p) were extended inside dK^ this set would 
represent the outer boundary of the largest set on M which, based solely on the 
geometry of (M, g, fc, J, p), must be contained within the black hole region of space 
time [9], [17]. Thus, although dK is, in general, not the set in which the spacelike 
slice intersects the event horizon, within the initial data set itself it provides the 
best possible representative of the black hole. In addition, if the censorship and 
energy assumptions above hold, any data inside dK will not affect the Cauchy 
development outside of the event horizon, and for this reason we exclude it from 
discussion. Finally, as shown by Hawking [8], dK will always consist of a finite 
union of topoligical spheres, and so there has been no loss of generality in the 
topological assumptions on M. 

The problem of the construction of black hole initial data with apparent hori- 
zon boundaries has received much attention in recent years. It would be impos- 
sible to cite every relevant work here, but the interested reader may see, for in- 
stance [1], [5], [6], In most of these works, the main tool for constructing 
black hole initial data with an apparent horizon boundary has been the conformal 
method, and considerable progress has been made. This method involves taking 
an initial "free" metric on R'^\iir and modifying it with a conformal factor to en- 
sure that the constraint equations are satisfied. The apparent horizon condition 
appears as a Neumann boundary condition on dK. However, with this method one 
has no control over the values of the conformal factor on the horizon; it is therefore 
impossible to specify the horizon geometry with this method. 

Since in the present work we want to prescribe the geometry of the apparent 
horizon, we use instead the parabolic equation for prescribing scalar curvature, 
for which the horizon geometry will be contained in the freely specifiable part of 
the metric. In [1], Bartnik does this for quasi-spherical metrics, which allows for 
the construction of apparent horizons that are intrinsically round. In the present 
work we relax this condition. Since this is a first attempt in that direction we 
are also going to make the simplifying assumptions that the horizon only has one 
component, the initial data set has moment of time symmetry fc = 0, and the local 
mass density p is compactly supported. This reduces the problem to the following: 
Given a compactly supported local mass density p > on M.^\BroiO), and a metric 
h on dBro{0) = S^, construct a metric g on R"^\_B,.(,(0), asymptotic to the standard 
flat metric 5ij, satisfying R{g) = IGnp, and .such that g\dBr.g{o) = ^• 

To attack this problem we make the standard spherical-polar identification of M 
with [ro, oo) x and construct g in the form 

(1) g = u^dr^ + r^7, 

where 7 vanishes on dr and g\s^ = r'^j. The parabolic scalar curvature equation, 
which we write in the next section, may be viewed as a second order parabolic 
partial differential equation for the radial component u of the metric as long as 
•f^^rdrjAB > —4. The component 7 appears in the coefficients along with the 
scalar curvature R. 
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To see how we can obtain the apparent horizon condition for a metric written in 
the form ([T]) note that the Schwarzschild metric cast into this form reads 



1 

1 ' 



where 7 is a fixed round metric on S^. Of course, at first glance it appears that 
the metric is singular at ro, but it is well known that an appropriate coordinate 
transformation shows that can actually be extended to the boundary as a smooth 
metric on a manifold with boundary. The blow-up in the normal component is, in 
fact, responsible for the surface r = tq being an apparent horizon for this metric. 
Indeed, this follows since the extrinsic curvature of the foliation spheres of a metric 
in the form ^ is 

r / 1 d"f 



and so for the Schwarzschild metric one has x = y/ry/r — ro7. Of course, in this 
case we actually know that Srg is a spacelike slice of the event horizon. 

Our approach to obtaining a horizon at r = tq will be to imitate the Schwarzschild 
case as closely as possible while striving for the greatest generality in the constructed 
metric. Namely, we shall construct metrics in the form 

^'^ , 2 2 

r 

by choosing 7 at the outset and using the parabolic scalar curvature equation to 
solve for v in order to get R{g) = Wirp > 0. Then the horizon geometry is just 
given by r^7 and is thus contained within the freely specifiable part of the metric. 
We may implement this procedure for a large class of 7. In particular, one has 

Main Theorem. Let h be a metric on §^ with positive Gauss curvature. Let e > 
and let 7 denote the fixed round metric on Let j(r) be a smooth family of 
metrics on [ro, 00) that satisfies, for some large r2 > ro + e, the following: 

7(r) = 7(ro), r G [ro,ro + £) 
7(r) = 7, r > r2 

djAB AB A 

^^^1 > -4 
or 

k(7) > 

Then for a compactly supported function p on M = [r^, 00) x satisfying Wnp < 
K, and which also vanishes on [ro,ro + e), there exists asymptotically flat time 
symmetric initial data {g,0) on M of mass density p such that g has the form 

2 2 
-dr + r 7, 



1 

where v is C°° , independent of r on [ro,ro + e), and bounded above and below by 
positive constants. The metric g has an apparent horizon at r = ro with intrinsic 
geometry (S^,/i); the horizon is, in fact, totally geodesic. 
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This theorem is an immediate corollary of a slight generalization to be proved at 
the end of the last section, in which we replace the conditions k > 0, Idnp < k with 
a more general, but more complicated, integral condition involving p, /f. A, k, K\ro- 

A few remarks are also in order concerning the other conditions on the free part 
of the metric 7. The first of these is obviously what allows us to specify the horizon 
geometry. The second is a technical condition, but only affects the generality of 7 
on the region [ro, rg + e], which can be made arbitrarily small; we shall refer to this 
region [ro, + e] x as the collar region. The third is a technical restriction that 
can be removed and replaced with an assumption on the decay of 7 to 7. However, 
doing so would require a detailed analysis of the asymptotic behavior, and this 
would detract from the purpose of the current work. The next to last condition 
serves the technical purpose that it causes the parabolic scalar curvature equation 
to be uniformly parabolic for r strictly greater than ro, but it also implies that the 
mean curvature of the foliation will be positive since the mean curvature of these 
spheres can, in general, be calculated to be 



This fact ensures that outside of r = ro there are no other compact surfaces satis- 
fying 9^ < 0, and so r = ro is an apparent horizon. 

The outline of the paper is as follows: In the next section we write the para- 
bolic scalar curvature equation and present some necessary results concerning this 
equation that were obtained in previous works. In the last section we prove the 
slight extension of our main theorem in two steps. In the first step we construct the 
metric on the collar region by solving an elliptic equation for v that we obtain from 
the parabolic scalar curvature equation by separation of variables; the conditions 
on the collar region are precisely what allows us to do this. This step is used to 
overcome the fact that the parabolicity of the parabolic scalar curvature equation 
must break down on the horizon. In the second step, which is contained in the 
proof of TheoremlH the metric exterior to the collar region is constructed using the 
results from the first section. 



The proof of the main theoremnis essentially a problem in prescribed scalar 
curvature. That is, given the "free" part of the metric 7, and the function p, we 
want to choose the remaining component u such that the metric g = u^dr"^ -\- r^7 
has scalar curvature R = IGnp. To do this, we use the parabolic scalar curvature 
equation, which relates u to i?, 7: 




2. The parabolic scalar curvature equation 




where Kj^ is the Gauss curvature of 7(r) and 
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A, B are used to denote components with respect to local coordinates (6*^, 6"^) on S^. 
For derivations of this equation see [1] , [15] , [13j . The equation has also been used 
successfully in [T3] [15]. For a discussion of the use of this equation to construct 
non-time-symmetric initial data see [2] , [2] , [H] ■ 

Note that \i H > Q and u is bounded above and below by positive constants 
then this equation is uniformly parabolic. As already noted, this is equivalent to 
the positivity of the mean curvature of the foliation spheres. Of course, this requires 
some pointwise a priori bounds on the solution u, but in the case that u is positive 
and bounded initially, these are easily obtained by the maximum principle under 
appropriate assumptions on the free data. The result of this is contained in the 
next lemma, whose proof can be found in [15j . [16j . and is a slight generalization 
of a result contained in [1 . The proof is repeated here since it yields simple but 
important bounds on the components of the constructed metric. To state this, given 
a function / on [ri,r2] x /*(r) = supp^gs f{r,p), f^{r) = infpgs2 f{r,p). 

Lemma 1. For r E I = [ri,r2\, let 7(r) he a family of metrics on §^ and let 
H, K, A, R be defined in terms of 7 as above. Define 



n 7,., \H\2 J J V7r, V J t 



ri 

Assume H > and K < 00. Suppose that u is a solution of equation ^ on [ri, r-2\ 
such that "initially" one has < u(ri, •) < 1/^/K. Then there exists a constant C 
depending on 7, R,ui, I such that C^^ < u < C. 

Proof. In order to prove the result, it is useful to introduce the auxiliary function 
w = to convert Equation ([2|) into an equation which is linear everywhere but 
the principal part: 

2 Au 2A 2 fr^R 

rdrW = — = ^ 

H u H H \ 2 

By using the maximum principle (see [16]) one can bound a solution of this equation 
with initial data wi from below by a solution of the ordinary differential equation 

dw^ ^2AV ( 2 fr^R 



dr \H ) \H\ 2 

with any initial data satisfying i(;:,(ri) < w{ri,-). For convenience put a = 2A/H,b = 
r'^R/2 — K so that the previous equation can be written rw' + a*w = —b*. Using 
the integrating factor exp — l)/tdt converts the equation into the form 

— rexp / \ dtw^ ] = —b cxp / dt, 



which can be immediately integrated to yield 

= - exp j ^ — -j_ — ^ dt i^riw^{ri) ~ J 6* exp J ^ — - — ^dtds 

The hypothesis of the lemma ensures that on the interval / the right hand side will 
be bounded below by a positive constant. This in turn gives the upper bound for 
u. To obtain the lower bound for u, we note that by using the maximum principle 
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again w can be bounded above by solutions of 

dw* /2A\ ^ ( 2 (r^R 



^ dr \H J ^"^ \H \ 2 

satisfying w*{ro) > w{ro,-). □ 

If, in addition, the coefficients of Equation ([2]) are C°° then the bounds on u 
just obtained guarantee existence of a C°° solution with the initial data ui. This 
is carried out in [T], |15| . The result is contained in the next theorem. 

Theorem 2. Suppose that the hypotheses of the previous lemma are satisfied for 
7, i? G C°^'i[ri,r2] x S^) andui e C°°{S'^). Then there exists a unique C°° solution 
u of Equation ^ on I = [ri,r2] with the initial data u(ri,-) = Ui. Furthermore 
C^^ <u < C, where C is the constant from the previous lemma. 

Note that in the previous theorem r2 can be as large as we would like so that, 
in fact, we get existence as r ^ oo as long as the free data continue to satisfy 
the hypotheses of the theorem on any interval. In our case we would actually like 
to obtain an asymptotically flat metric as a consequence of this. Hence, we need 
to assume appropriate asymptotic behavior on 7, i?, which in our case we have 
done, and derive from this the correct asymptotic behavior of the solution u. As 
in previous work for Equation ([2]) , this is done by proving the boundedness of the 
function m defined by 

o 2m 

u-^ = 1 . 

r 

The function m verifies the equation 

A 7; 

HrdrTu = — !- - {2A- H)m + r{A - B). 

Since we have made the assumptions that on a neighborhood of infinity 7 is identi- 
cally the round metric and i? vanishes, the equation simplifies dramatically. Indeed, 
we have H = 2,A—\ so that the equation for m becomes 

(3) 2rdrm=^. 

u 

Bartnik has made a detailed analysis of this equation in [IJ, in which the following 
is proved: 

Theorem 3. Suppose that m is defined and C°° on [ro,cxD) x §^ and satisfies 
Equation ([3|), where m,u are related as above. Then there is a constant toq and 
G §^ verifying (A + 2)ip ~ such that 

m = Too H h e{r), 

r — 2too 

where e is a C°° function on [ro,oo) x satisfying, 

\{rdryV^e\ < Qjr^Mogr, 
for all i, j > such that i + j < k. 



This theorem shows, in particular, that the constructed metric g will be asymp- 
totically flat in the sense mentioned in the introduction. Furthermore, the mass of 
the manifold is toq. 

Putting the previous results together, we get the result needed to construct the 
metric outside of the collar region. 
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Theorem 4. For r G [ri , cxd) let j{r) be a family of metrics on satisfying H > 
and K < oo, where K — if[ro,oo) defined as in Theorem\^ Furthermore, assume 
for some r2 > ri that 7 is the fixed round metric on and R=0 for r > r2 . Then 
if ui £ C°°(S^) satisfies < ui < \ j\fK, there exists a unique C°° asymptotically 
flat metric g in the form 

2j 2 I 2 
g — u dr + 7" 7, 

with u{ri, •) — ui, and whose scalar curvature satisfies R{g) = R- 

3. CONSTRUCTING THE DATA 

To construct the black hole initial data we construct the data on the collar 
region, and make sure this patches smoothly with data constructed via Theorem |4] 
for r > ri. At first glance, it seems problematic to use Equation ^ to solve for u 
on the collar region since at r = rg the equation must fail to be parabolic if this 
surface is to be a horizon. This problem is resolved by assuming that 7 is fixed 
on this region. Assuming in addition that v is also fixed in r on the collar region 
yields an elliptic equation for v on S^. Indeed, the fact that 7 does not vary in r 
immediately gives that II = 2, A = 1 so that Equation ([2]) becomes 

2r— = u A~fU + u — Ku . 
or 

Substituting u — v {1 — r^/r) ^ yields the following equation for v. 
(4) A^-y - Kv +- = 0. 

V 

One has 

Theorem 5. Assume k > and let = infs2 k, k* = supg2 k. Then Equation (|4]) 
has a positive solution v G C°° satisfying l/^f^F < v < I/^/k^- Within the class of 
functions there are exactly two solutions ztv. 

Proof. We are going to produce the solution by the method of sub-solutions. Before 
starting the process, we note that sub(super)-solutions for Equation ([4]) in the sense 
of 

A^v - Kv + -> (<)0 

V 

are also sub-solutions and super-solutions in the sense that any super-solution 
bounds any sub-solution from above. As usual, this is just a consequence of the 
maximum principle, but we should carefully check this anyway: Assuming that 
v^,,v* are sub and super solutions of Equation Q in the sense of the previous 
inequality, we consider the equation for the difference Sv* = v* — 

ASv* - {k+^)Sv* < 0. 

To see that Sv* > 0, suppose instead that Sv* < somewhere, and consider the 
set = {p e : 6v*{p) < 0}. Since the term in parentheses is positive, and since 
we have assumed fi is nonempty, by the maximum principle one has infn Sv* = 
infan Sv* , which is a contradiction since Sv* — on dil. 
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To now begin the method of sub-solutions, we are going to take as our starting 
sub-solution = For future reference, note that v* = I-^kT is a super- 

solution. Choose A > and define Vi as the solution of 

Avi — (A -f k)vi — — Au* . 

Then an application of the maximum principle just as in the first paragraph shows 
that vi > since the difference 6vi — vi — satisfies 

ASvi - {\ + k)Svi < 0. 

We now define Vi inductively as the solutions of 

(5) Avi ~ {X + K)vi ^ -Xvi-i —. 

Vi-l 

The fact that Vi > Vi-i now follows from induction and another application of the 
maximum principle since Vi-i > Vi-2 yields for the difference Svi = Vi — Vi-i the 
inequality 

f 



ASvt - (A + k)Svi = - a 6vi-i < 0. 

V V^-lVi^2J 

Thus, we now have an increasing sequence of functions Vi, which are, in fact 
sub-solutions since for the right hand side of ^ we have 

-Xvi-i — = -\v, h 1 A ) {vi - Vi^i) > -Xv, . 

Vi-l Vi \ V^Vi-l ) Vi 

Hence, we now have a pointwise non-decreasing sequence of sub-solutions Vi that 
satisfies v^, < Vi < v* . Thus Vi converges to some function v pointwise. In addition, 
the bounds on Vi imply that the right hand side of the Equations ([5]) are uniformly 
bounded in i'' for any q. Whence, by applying elliptic regularity theory we see that 
the Vi are also uniformly bounded in VF^'*. Applying elliptic regularity iteratively, 
we may bootstrap to see that, in fact, the Vi are uniformly bounded in W'^''^ for any 
k,q so that the Sobolev embedding theorem shows that they are uniformly in C'^ 
for any k. The Ascoli-Arzela theorem then shows that a subsequence converges in 
C*^ to a function that can be none other than v. Taking the limit in Equation ^ 
then shows that v must be a solution of Equation ([4]) . 

Concerning the uniqueness, first note that if a solution is to be of class C^, then 
it cannot change signs, and so we may restrict to the class of positive functions. 
Given two positive solutions vi,V2 one can see that vi — V2 by, again, applying the 
maximum principle to the equation for the difference Sv = V2 ^ vi 

A6v- ( K + ^— ) Sv = 0. 
V V1V2J 

□ 

We may now prove our final theorem, of which the main theorem is an immediate 
corollary. 

Theorem 6. Let ^ be a family of metrics as in the hypothesis of the main theorem 
with the exception that we do not assume ^(7) > outside of the collar region. 
Instead, assume that 7, p are such that, with R = Wnp and K defined as in Theo- 
rem^ one has 

K < rlK{h) (1-- 
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Then there exists asymptotically fiat time symmetric initial data {g, 0) on M of 
mass density p such that g has the form 

""^ , 2 2 

r 

where v is C°° and hounded above and below by positive constants. 

Proof. From the previous theorem we have the metric in the above form on the 
collar region. Taking ui — w/-^/! — ro/ri in Theorem |4] gives the metric for r > ri 
since < ui < 1/ We only have to check that the constructed metric is smooth 
across r = ri. This is easily done by applying Theorem |4] again for r G [vi — S,oo) 
for very small 6. Indeed, the new metric is certainly C°° at r = ri and agrees with 
the previous metric on [ri — d^oo) by the uniqueness part of that theorem. □ 
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